Invertible and non-invertible alloy Ising problems 



C. Wolverton and Alex Zunger 
National Renewable Energy Laboratory, Golden, CO 80401 

B. Schdnfeld 

Institut fur Angewandte Physik, Eidgenossische Technische Hochschule Zurich, CH-8093 Zurich, Switzerland 

(February 6, 2008) 



Physical properties of alloys are compared as computed from "direct" and "inverse" procedures. The 
direct procedure involves Monte Carlo simulations of a set of local density approximation (LDA)- 
derived pair and multibody interactions {fj}, generating short-range order (SRO), ground states, 
order-disorder transition temperatures, and structural energy differences. The inverse procedure 
involves "inverting" the SRO generated from {v/} via inverse- Monte-Carlo to obtain a set of pair 
only interactions {£'/}. The physical properties generated from are then compared with those 
from {vf}. We find that (i) inversion of the SRO is possible (even when {vf} contains multibody 
interactions but {£/} does not) but, (ii) the resulting interactions {£7} agree with the input inter- 
actions {vf} only when the problem is dominated by pair interactions. Otherwise, {£7} are very 
different from {vf}. (iii) The same SRO pattern can be produced by drastically different sets {17}. 
Thus, the effective interactions deduced from inverting SRO are not unique, (iv) Inverting SRO 
always misses configuration-independent (but composition-dependent) energies such as the volume 
deformation energy G(x); consequently, the ensuing {£7} cannot be used to describe formation 
enthalpies or two-phase regions of the phase diagram, which depend on G(x). 



PACS numbers: 7I.I0.+X, 64.60. Cn, 64.70.Kb 



The physical properties of A\- X B X alloys are usu- 
ally analyzed and interpreted via "cluster expansion" 
modelsrlM Each of the N sites of an alloy i=l,...,N is 
labeled by a spin variable Si = -1 or +1 if site i is oc- 
cupied by an A or B atom, respectively. The set of spin 
variables {Si} defines a configuration a. The energy of 
any of the 2 N possible configurations is then written as 
a sum over clusters of points {i; ij; ijk; ...}:□ 



E(a, V) = J (V) + J2 UV)^ + J2 J ii{V)SiSj 

i j<i 

+ ^ Jijk{V)SiSjSk + •■• 
k<j<i 



(1) 



where V is the volume, the J's are interaction energies 
and the first sum is over all sites, the second over all 
pairs, the third over all triplets, etc. We refer to these 
elementary clusters as "figures" /. 

If the set of interactions {Jf(V)} is known for a given 
alloy system, one may apply standard methods of lat- 
tice statistical mechanics (e.g., mean field, cluster vari- 
ation, or Monte Carlo methods) to the expansion and 
compute ground state structures or finite-temperature 
thermodynamic properties. Recent examples include the 
calculation of temperature-composition phase diagrams 
and ground state structures of transition metaltla'B and 
semiconductor^ alloys, mixing enthalpies of disordered, 
partially ordered, and off-stoichiometric alloysjim and 
short-range order (SRO) of solid solutions BO We refer 
to this approach as the "direct approach". 

Conversely, another common tradition involves the 
"inverse approach": A measured thermodynamic prop- 



erty such as the set of SRO parameters a(n) (the 
atom-atom pair correlation for the nth atomic shell) is 
used in an inverse statistical approach (e.g., the inverse 
Monte Carlcu/IMC) methodQ) to deduce a set of effective 
interactions .u These interactions are subsequently used 
in a cluster expansion [Eq. ([!])] to predict thermody- 
namic properties other than the SRO. In this paper, we 
explore the extent to which the inverse approach may be 
used to predict alloy properties by applying it to a well- 
characterized a(n) obtained through a direct procedure. 

In the following it is convenient to introduce the excess 
energy AE(a, V) of configuration a defined as the energy 
of this configuration at volume V, relative to the energies 
Ea(Va) and Eb(Vb) of equivalent amounts of solid A and 
B, at their respective equilibrium volumes Va and Vb'- 

AE{a, V) = E(a, V) - [(1 - x)E A {V A ) + xE B {V B )]. (2) 

If the equilibrium volume V(a) depends primarily on the 
composition x and only weakly on the configuration c, 
then the variables a and x can be rigorously separated 
in Eq. (0) givingp 



AE(o-,V) = G(x) 



f 



(3) 



Here, the first term G(x) describes the elastic energy nec- 
essary to deform the constituents from their equilibrium 
volume to the volume V(x) of a. The second term of 
Eq. (||) describes the spin flip excess energy of forming 
a from A + B already prepared at the volume V. The 
correlation function 11/ is defined as a product of the 
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TABLE I. The values of the input interaction energies Vf and the interaction energies iif reconstructed via IMC simulations 
of the SRO computed from Vf (meV/atom). Designation of the "figures" / follows the notation of Table IV of Ref. [| For 
Set 2, the multibody interactions used in the direct set DfUf, but not in the inverse set DjVf are (in meV/atom): Js=— 96.1, 
K 3 =44.5, 1/3=64.5, M$=— 41.1, Qs~ — 81.3, and _K" 4 =139.1. Structural energy differences, ordering energies 8E OTd (a) (the 
energy difference between a and a random alloy at the same composition), and the random alloy mixing energy at a;=l/4 are 
shown (meV/atom), as are transition temperatures (K). "NA" means not applicable. 

Set 1 Set 2 







Direct 


Inverse 


Direct 


Inverse 


Clusters 


Designation 


DfUf 


D f i> f 


DfUf 


Df i/f 


Empty 


Jo 


-233.2 


NA 


-233.2 


NA 


Point 


Ji 


252.9 


NA 


252.9 


NA 


Pairs 


■h 


152.0 


157.2 


152.0 


690.0 




K-2 


-20.0 


-21.0 


-20.0 


17.6 




L-2 


58.9 


60.0 


58.9 


-19.2 




M 2 


33.5 


33.3 


33.5 


103.2 




N 2 






0.0 


-3.6 




Oi 






0.0 


-0.4 




P2 






0.0 


13.2 




5E{Ll 2y D0 22 ) 


-4.0 


-4.8 


+103.3 


+76.2 




SE old {Ll 2 ) 


-42.6 


-45.1 


-41.9 


-79.8 




Ground State 


Ha 


Ll 2 


D0 22 


-DO22 




T c 


630 


680 


1850 


1900 




G{x = 1/4) 


112.1 


NA 


112.1 


NA 




A# mix (l/4) 


-56.2 


-172.1 


-227.7 


-600.6 



variables Si over all sites of the figure / with the overbar 
denoting an average over the Df symmetry equivalent 
figures per lattice site. Equation (|3|) is similar to Eq. 
(|l|), but here the effective interaction energies {y/} are 
volume-independent pure spin flip energies. 

We will examine the invertibility of the inverse ap- 
proach by performing a "controlled experiment" : As in- 
put, we use two "exact" sets of interactions, {v;}. Equa- 
tion (|l|) is then used along with these {vf} in direct 
Monte Carlo (MC) simulations to obtain the "exact" 
quantities such as SRO parameters a(n), structural en- 
ergy differences 6E(o~,o~') between configurations er and 
a', transition temperaures T c , and the mixing energy of 
the random alloy AiJ m j X . We then contrast the results 
of this "direct procedure" with those of the "inverse pro- 
cedure", in which the set {a(n)} (obtained in the direct 
procedure from the known, exact {ff}) is used as input 
to deduce the interactions {£/} by IMC simulations from 
which we then obtain a(n), 6E(a,a'), T c , and AH mlx . 

We use as input two sets of interaction energies {vf} 
(see Table |). As an illustration of physically realistic 
interactions, wejise one set of interactions that was re- 
cently extracted!! from T — first-principles calculations 
of formation energies of ordered fcc-based Nii_ x V x com- 
pounds and reproduces reasonably well many of the mea- 
sured physical properties. This "realistic set" , which we 
call Set 2, contains pair interactions up to fourth neigh- 
bors, as well as three and four-body interactions. Set 1 is 
identical to Set 2, except that we have set equal to zero 
all multibody interactions. 

We first contrast the directly calculated alloy proper- 
ties using Sets 1 and 2 in MC simulations. For the direct 
MC calculations, a system size of 4096 atoms was used 



with periodic boundary conditions, 1200 Monte Carlo 
steps (MCS) were used for equilibration, and averages 
were typically taken over 1800 MCS. Temperatures of 
T=850 K and 2300 K were used for the SRO calcula- 
tions for Sets 1 and 2. Fig. [I] shows the SRO a(k) 
calculated directly from {f/} for Sets 1 and 2 at com- 
position A3B. Only multibody interactions contribute 
to the difference between the SRO of Sets 1 and 2, and 
this difference is dramatic: a(k) of Set 1 shows peaks at 
the X-points (100) whereas a(k) of Set 2 shows peaks at 
the W-points (1 i0) as seen experimentally in N13V.E2I Ta- 
ble H also shows that multibody interactions change the 
ground state structure from LI2 to the observeat3 DQ22 
structure and that the energy difference between these 
two structures, 5E(L1 2 , DO22), changes from -4 to +103 
meV/atom and T c changes from 630 to 1850 K upon in- 
clusion of multibody interactions. Also, note from Table 
I that G(x = 1/4) is a significant fraction of the random 
alloy mixing enthalpy, A_ff m j x (a; = 1/4). Thus, from the 
directly calculated values it is clear that both multibody 
interactions and the elastic energy G{x) are physically 
very important in this alloy system. 

Using the directly calculated a(k), we now apply IMC 
to recover the interactions energies. Following the tradi- 
tion among practitioners of the IMC method, only pair 
interactions were retained in the energy expression of 
IMC. First, configurations were produced which repro- 
duced the input values of 35 shells of a(n). System sizes 
of 262,144 and 216,000 sites were used for Sets 1 and 2. 
IMC simulations were performed on three crystals com- 
patible with the sets of a(n), and averages were taken 
over these three crystals. Tests were performed of the 
convergence of the inverse procedure with respect to the 
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FIG. 1. Comparison of the (001) planes of a(k) and a(k). 



number of pairs included: Many different sets of pair in- 
teractions were used (between 4 and 20 shells) , and from 
these calculations, the number of pairs needed to ade- 
quately reproduce a(n) were determined. 

Short-Range Order: Fig. [j] compares the recalcu- 
lated a(k) (computed from vf) and a(k) (computed 
from Vf). Both sets of SRO are well reproduced by 
the IMC procedure. The accurate inversion of SRO [i.e., 
{a(n)} {a(n)}] has been demonstrated previously by 
many authors (see, e.g., Refs. fj], ||). However, in these 
previous studies, measured SRO was used as input to the 
IMC, and thus the interactions which produced the input 
SRO were not known. We have shown that even when 
multibody interactions are used to produce a(n) (as in 
Set 2), the IMC procedure (using pair interactions only) 
reproduces this SRO quite well. 

Effective Interactions: We compare the values of Vf 
vs. Vf in Table Q. For Set 1, the IMC algorithm closely 
reproduces the input set of pair interactions; The stan- 
dard deviation between {vf} and {£/} for Set 1 is 2.7 
meV/atom. For Set 2, the direct and inverse sets of inter- 
actions differ dramatically in several respects: (i) Three 
and four body interactions are present in {vf}, but are 
specifically excluded from the pair-only IMC calculation 
of {£/}. (ii) There are huge differences in pair interac- 
tions (the standard deviation of the first four pairs is 275 
meV/atom). (hi) Seven pair interactions were required 
in the IMC to reproduce a(n), whereas only four pair 
interactions were present in the direct set. (A similar in- 
crease in range was reported in Ref. |l^.) However, even 
though there are enormous differences between Vf and Vf 
of Set 2, they both produce nearly identical SRO patterns 
(see Fig. [j]). This surprising fact indicates that even 
the pair interactions are not determined uniquely from 
a SRO pattern. [This non-uniqueness was also found by 
Schweika and Carlsson (Ref. [l^; Fig. 3a), who in contrast 
to the present work, used a high-temperature expansion 
(whereas we use IMC) to invert SRO generated by pair 



and multibody interactions.] We assert that due to the 
non-uniqucness of pair interactions derived from IMC, 
they cannot generally be compared with other sets of pair 
interactions. When multibody interactions are physically 
important, the non-uniqueness of these sets make such 
comparisons meaningless. For example, Schweika and 
CarlssorJld found that inversion of SRO produced interac- 
tions Vf which were temperature- dependent even though 
the input set {vf} was not. Clearly, this temperature- 
dependence is not due to physical effects (e.g., vibrational 
or electronic excitation effects), but rather due to the fact 
that a pair-only inverse scheme does not recover informa- 
tion on the multibody interactions {vf}. 

Structural energies, ground states, and transition tem- 
peratures computed from {vf} are compared with those 
computed from {£/} in Table [jj 

Mixing Energies: Table [| also shows values of the mix- 
ing enthalpy of the random alloy AH m i K (x = 1/4). For 
Set 1, where {£/} w the direct and inverse val- 

ues of AiJ m j X differ by more than 100 meV/atom, even 
in cases (such as Set 1) dominated by pair interactions. 
In Set 2, the comparison of A_ff m i X is even worse (direct 
and inverse values differ by more than 350 meV/atom). 
Deducing values of AiJ m ; x is clearly not reliable in the 
inverse procedure. 

We have seen that while that total energy E(a, V) de- 
fines the complete set of interaction energies, inversion 
of quantities (e.g., SRO) other than E(a, V) may lead to 
a loss of information. We now use Eq. (||) to distinguish 
different classes of alloy properties and discuss which are 
invertible: 

(a) Physical properties that depend on both G(x) and 
on the spin flip energies {vf} include any quantity which 
involves the energetics of two or more concentrations and 
hence, two or more volumes. [Note that G(x) depends on 
x, but not on the particular atomic arrangement ("con- 
figuration" ) a.] Examples include the formation energy 
of a structure [which involves V, Va, and Vb, c.f. Eq. 
(0)], the mixing energy A_ff m j X of the random alloy, and 
two-phase equilibria in a composition-temperature phase 
diagram. Since "type- (a)" properties such as the set 
{AE(o-, V)} contain complete information on both G(x) 
and on all {vf}, given the measured or ab-initio calcu- 
lated energies {AE(a : V)}, it is possible to invert Eq. (|3|) 
and in principle extract the "exact" G{x) and {V/}, as 
demonstrated in Ref. |9|. Thus, inversion of "type-(a)" 
properties involves no loss of information. 

(b) Physical properties that do not depend on G(x) 
include energy differences of iso-compositional configu- 
rations cr and a', 8E(o~,o~'). The order-disorder transi- 
tion temperature T c at stoichiometric composition also 
falls into this class since it involves the energy differ- 
ence between the disordered high-temperature phase and 
the partially ordered low-temperature phase, both at the 
same x. Another physical property which does not de- 
pend on G{x) is the atomic SRO. a(n) involves a com- 
petition between energies of a random and a short-range 
ordered structure, both at the same volume V(x)\ there- 
fore, a(n) (even if determined for several compositions) 
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contains no information about G(x). Therefore, inver- 
sion of a "type-(b)" property, such as SRO, cannot pro- 
vide any information on G{x), even if the SRO covers 
a range of compositions. Consequently, the interactions 
{£/} extracted from such a procedure do not allow calcu- 
lation of "type-(a)" properties, such as formation energies 
[Eq. (Q)], mixing energies, or the phase-coexistence re- 
gions of the phase diagram. This point is highlighted by 
recent studiesQ on Nii-^Au^: This is a phase-separating 
system; however, the SRO is of ordering type. Inverting 
the SRO will thus inevitably produce ordering-type {£/}, 
which are useless for predicting the miscibility gap phase 
diagram-|Or the correct AiJ m j X > 0. These conflicts are 
resolvedB by using G{x) in the Ising-like expansion. 

We conclude that: (i) The IMC algorithm provides a 
set of pair interactions which accurately reproduces the 
input SRO whether or not multibody interactions are used 
to generate this input. When only pair interactions are 
involved, the inverse procedure can even provide accurate 
values of structural energy differences, ordering energies, 
and energies of SRO; However, (ii) when multibody inter- 
actions are physically important, even the pair interac- 
tions are incorrectly determined by the inversion of SRO. 
The structural or ordering energies deduced from the in- 
verse procedure can thus contain substantial errors, (iii) 
Finding a set of interactions which reproduces a given 
set of SRO is found to be a non- unique process: dra- 
matically different sets of interactions (one set with pairs 
only, one set with pairs and multibodies) may still pro- 
duce quantitatively the same SRO. Thus, comparing sets 
of interactions from IMC with other sets of interactions 
may be unwarranted. However, comparing a theoretical 
SRO pattern to a measured one is a sound procedure, 
(iv) Inverting the SRO always removes information on 
energy terms that are SRO-independent, e.g., G(x). This 
loss prevents, in principle, the interactions deduced from 
SRO from being applied to predict phase-coexistence re- 
gions of the phase diagram or AiJ m ; x . 
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